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I. INTRODUCTION 

About a decade ago, in a remarkable series of papers, 
Tan showed that high-momentum correlations, in few- 
and many-body quantum systems with short-range in¬ 
teractions, are governed by a single number: the con¬ 
tact C P]. In particular, he found analytically that 
the high-momentum tail of the momentum distribution 
obeyed the law n(k) ~ C/k 4 , and that C also governed 
the behavior of the energy upon small adiabatic changes 
in the coupling strength (as driven, e.g., by an exter¬ 
nal magnetic field in ultracold atom experiments M). 
These developments took place more or less in parallel 
with the work of others (see, e.g., Ref. a ), and they 
were eventually understood in terms of the operator- 
product expansion of quantum field theory by Braaten 
and colleagues [5j . Other works followed that considered 
the high-momentum behavior of more complex quantities 
(such as the stress energy tensor), which led to the ap¬ 
pearance of the contact in sum rules involving the shear 
and bulk viscosities, the superfluid density, and other re¬ 
sponse functions EH]. 

While those advances were by and large analytic, the 
calculation of C itself in few- or many-body systems 
typically requires solving those problems numerically, 
e.g., via quantum Monte Carlo methods. The large- 
momentum and adiabatic relations mentioned above (see 
also Refs. mm), however, provide possible avenues for 
the computation of C, which will in general depend on 
all the dynamic and thermodynamic variables of the sys¬ 
tem. Calculations of C using nonperturbative methods 
(Bethe ansatz in one dimension, lattice and diffusion 
Monte Carlo methods in two and three dimensions) ap¬ 
peared in Ref. m in three dimensions, in Ref. [12] in two 
dimensions, and in Ref. m in one dimension. Excellent 
reviews can be found in Refs. mm- 

In this work, we present our calculations of the 
ground- and first-excited-state energy and Tan’s contact 
for multiparticle nonrelativistic fermion systems in a one¬ 
dimensional (ID) box with periodic boundary conditions 
(i.e., a ring). We cover the range from few- to many-body 
systems and couplings from noninteracting to strongly 


coupled. On the experimental side, this problem has been 
explored in Ref. d, although that realization differs 
from the present system as it features an external har¬ 
monic trapping potential. On the theory side, the prob¬ 
lem can be solved using the Bethe ansatz technique OBI, 
but we have chosen to approach it using lattice Monte 
Carlo methods (i.e., lattice field theory formulations pow¬ 
ered computationally by lattice-QCD techniques). The 
latter are applicable in situations where the former is 
not, notably in higher dimensions and in the presence of 
an external potential. Furthermore, lattice methods are 
an area in which remarkable strides are currently being 
made across physics: from materials science m to nu¬ 
clear physics HD], and of course lattice QCD [21]. There¬ 
fore, understanding the advantages, disadvantages, and 
potential of these methods is both interesting and timely. 

II. HAMILTONIAN, SCALES, AND 
DIMENSIONLESS PARAMETERS 

As anticipated above, we focus on a ID system of non¬ 
relativistic fermions with a zero-range interaction (also 
known as the Gaudin-Yang model [22]) and periodic 
boundary conditions, such that the Hamiltonian is given 

by 

^ = -£X V ?-Z 9S( Xi - Xj ), (1) 

i i<j 

where the sums are over all particles. We restrict our¬ 
selves to a two-species unpolarized system, but higher 
degeneracies can be studied with the same methods. We 
employed the same technique as in Refs. [23] adapted 
to one dimensions. Following that approach, we placed 
our system in a Euclidean spacetime lattice of extent 
N x x N t and used a Trotter-Suzuki decomposition of the 
Boltzmann weight followed by a Hubbard-Stratonovich 
transformation jM]. The path integral is evaluated us¬ 
ing Metropolis-based Monte Carlo methods (see, e.g., 
Ref. [25]), specifically in the form of the hybrid Monte 
Carlo algorithm. Further details on the method em¬ 
ployed here are explained below. 
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In the following, we use units such that K = m = 
ks = 1, where m is the mass of the fermions. The 
physical input parameters are the total particle number 
N = Ah + iVp the size of the box L = N x £ (where t = 1 
to set the length and momentum scales), and the (at¬ 
tractive) coupling strength g\ only the last two of which 
are dimensionful. We use all of these to form one dimen¬ 
sionless intensive quantity: the dimensionless coupling 7 
given by 


7 = g/n, ( 2 ) 

where n = N/L is the particle-number density, as is 
common in other ID ground-state studies (see, e.g., 
Refs. [21127]). The extent of the temporal direction is 
/? = tN t , which we vary in order to extrapolate to the 
large-/3e F limit, where e F = k F /(2m) and k F = nir/2; 
we further elaborate on the relevant scales below. Note 
that, in one dimension, fermions with a contact interac¬ 
tion arc ultraviolet-finite, and therefore the bare coupling 
has a well-defined physical meaning, namely, g = 2 /a 0 , 
where a 0 is the scattering length for the symmetric chan¬ 
nel (see, e.g., Ref. [28]h 


III. MANY-BODY METHOD 


The (unnormalized) ground state of a many-body 
quantum system with Hamiltonian H can be obtained 
as the large-/3 limit of 

W)p = ex P (-/?#) IV'o). (3) 


as long as the “guess” state | 4 ) has a nonvanishing pro¬ 
jection onto the true ground state | ip). Thus, one may 
write the ground-state expectation value of an operator 
O as 


(6) = lim Ojj , 

p—± OO 


where we have defined 

Q _ (i’ 0 \U(/3,P/2)dU(/3/2,0)\^ 0 ) 
0 <4|t/(/3,0)|4> 

and the imaginary-time evolution operator 


U(n,T a ) = exp 


-{Tb - T a )H 


(4) 

(5) 

( 6 ) 


In this work we will use this simple formalism and take 
| 4 ) to be a Slater determinant of single-particle orbitals 
{4> k } given by plane waves, where k = 1 , 2 ,..., TV*, with 
Atj. = = N/2 being the number of fermions of each 

flavor. 

We approximate the operator U using a Suzuki-Trotter 
decomposition 


— rf/2 e —W e — tT/2^^ (7) 


where r is our imaginary-time discretization parameters, 
and we have split the Hamiltonian into the one-body ki¬ 
netic energy operator T and the two-body, zero-range 
interaction V . At each time step t , we implement an 
auxiliary field transformation writing (generically) 

e- T * = J Va(x)e- T9 t.«e~ T9 ^ ( 8 ) 


where the V s a are one-body operators that depend on the 
Hubbard-Stratonovich field <r(x) and f Va(x) is a sum 
over all possible configurations of a at the specific time 
slice t. After collecting the factors corresponding to each 
time step, we identify the zero-temperature partition sum 

(V’ol U(/3,0) \i/j 0 ) = J Va(x,t) P[a], (9) 


where now the path integral is over a spacetime varying 
field cr(x,t), and we defined the Monte Carlo probability 

P[a] = (%\ U a {fi, 0) |4) = det 2 [M„m , (10) 


where 


U*(T a +T,T a ) = e rT/2 e rV t.-e TV ^e rT/2 , 


( 11 ) 


and the matrix M a ([3) is the single-particle representa¬ 
tion of the product operator U„{&,$), restricted to the 
sub-space of the Hilbert space spanned by the orbitals 
{<t>k}, i-e., 


[M a ((3)] ab = (a\UAP,m. ( 12 ) 

The square of the determinant in Eq. ( |10| ) results from 
the fact that we are considering two (distinguishable but 
otherwise identical) fermion species. Sampling the auxil¬ 
iary field according to P[u], one evaluates the expectation 
value of observables using 

Op = J T>aP[a]0[a), (13) 

as a function of the imaginary time /3, where 

^ Wo\U'rW,W)dU'rWM\'l’o) (M) 

( 41 ^ 03 , 0 ) 14 ) 

This is followed by an extrapolation to large /3. 

Using this technique, we compute the ground-state en¬ 
ergy Eq and Tan’s contact C. To determine the ground- 
state energy, one of the simplest ways to arrive at the 
observable of interest is actually to differentiate InZ, us¬ 
ing Eq. with respect to r. This is equivalent to us¬ 
ing Wick’s theorem [25], and it generates a simple and 
normalized expression for a stochastic estimator of the 
energy, 



E* 


din Z 

<9/3 



Va P[a] tr 


M~\P) 


dM a (i 3)1 

dr 


U[T a + T,T a ) = e 


(15) 
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In the above expression, upon performing the last differ¬ 
entiation, it is easy to identify the contribution to Eg 
coining from the kinetic and interaction energies: us¬ 
ing Eq. (12 1 , the r derivative will act on all the factors 
present in U and bring down kinetic- or potential-energy 
terms from the exponent. 

Using the estimator for the potential energy, we can 
calculate Tan’s contact C, which is defined by 


C = 2 -^ - = —g(V), (16) 

da 0 

where we have used the Feynman-Hellman theorem and 
the fact that V is just a contact interaction, which sat¬ 
isfies g dV/dg = V. The above definition of C holds for 
the ground state as well as for excited states. 

In the limit of large /3ep, the difference between 
the finite-imaginary-time expectation value Ep and the 
ground-state energy decays exponentially at a rate deter¬ 
mined by the difference between that minimum energy 
and the energy of the first excited state E\: 

Ep = E 0 + Ae~^ El ~ Ea) + O f e -^ E2 ~ E °A . (17) 


From this decay, we determine the energy of the first 
excited state as a function of the coupling 7 , and proceed 
to compute the excited state contact C\ via 


C i = 2 


d E 1 
da 0 



(18) 


To our knowledge, this is the first determination of the 
contact in an excited state, other than calculations at 
finite temperature. In the next section we present our 
results for these quantities. 


IV. ANALYSIS AND RESULTS 
A. Ground- and excited-state energies 

We performed calculations for N = 4 — 12 fermions in 
lattices of spatial sizes up to N x = 40 —100 and couplings 
7 = 0, 0.2, ..., 4.0. We used r = 0.025, which resulted 
in temporal lattice sizes varying between 60 and up to 
10 4 points. For each set of parameter values, we took ap¬ 
proximately 5000 decorrelated samples of the auxiliary 
field a (which yields statistical uncertainties on the or¬ 
der of 1 - 2 %) and evaluated the total-energy estimator of 
Eq-©- Extrapolating to the large-/3e f limit (i.e., large 
N t limit), we obtained the ground-state values. In each 
case, studying the exponential decay of this estimator for 
increasing imaginary time allowed us to obtain the energy 
of the first excited state. In Fig. [l] (top panel) we show 
typical Monte Carlo results for the total energy estima¬ 
tor, with the corresponding fits used to obtain ground- 
and excited-state estimates. The error bars represent the 
statistical uncertainty of the Monte Carlo calculation. 



P e F 



Figure 1. (color online) Top: Estimate of the energy of 
12 unpolarized fermions in a N x = 40 lattice with periodic 
boundary conditions, in units of the noninteracting energy 
E fg = |£ f N = 7t 2 A 3 /(24Z/ 2 ), as a function of imaginary 
time fie F . From top to bottom, the lines correspond to cou¬ 
plings 7 = 0, 0.2, ..., 4.0. The data are shown with statistical 
error bars. The continuous lines show the least-squares fits 
for exponential decay to the ground state. Bottom: Estimate 
of the ground- and excited-state energies of 12 unpolarized 
fermions in a periodic ID box. The statistical uncertainties 
are smaller than the size of the symbols. 


The energy of the ground- and first-excited states, as 
a function of the coupling 7 , are shown in Fig. [l] (bottom 
panel) for a representative system. The error bars for 
the ground-state energy are taken from the exponential- 
decay fit used to extract the asymptote at large imagi¬ 
nary time (top of Fig. [l]). In the same figure we show the 
first-excited state energy. For the latter, the data were 
obtained by studying the exponential decay of Fig. [T] for 
all values of 7 and incorporating, using a power-law fit as 
a function of 7 , the exact value at 7 = 0. The associated 
uncertainties represent 95% confidence intervals. 

The above lattice results are typical. In order to obtain 
the physical values in the continuum limit, we extrapo¬ 
lated to N x —> 00 . The extrapolated data are shown in 
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Figure 2. (Color online) Top: Estimate of Tan’s contact in 
units of Lkp for N = 12 and N x = 40 for various couplings 7 
showing the extrapolation to large /3e f . From bottom to top, 
the lines correspond to couplings 7 = 0, 0.2, ..., 4.0. Bottom: 
Estimates of Tan’s contact for 12 unpolarized fermions in a 
N x = 40 lattice with periodic boundary conditions. The blue 
line shows the ground-state result, while the red line shows 
the excited-state value. 


Table [T] and further details on the extrapolation proce¬ 
dure are given in Sec. |VJ 


interaction energy, and the uncertainties shown were ob¬ 
tained in the same fashion. The contact associated with 
the excited state was computed via Eq. (18). As with 
the excited-state energy, the uncertainties in C 1 were ob¬ 
tained by performing a fit to our whole data set as a func¬ 
tion of 7 that accounts for the C 1 = 0 value at 7 = 0. 
Our full results for the contact, as a function of 7 and 
N, extrapolated to the infinite-volume limit, are shown 
in Table [HI 


C. Approach to the thermodynamic limit 

As can be appreciated in Tables |T| and |TT} the variation 
in E/E fg and C/{Lkp) as a function of particle number 
is relatively small. This behavior is indicative of a rather 
fast approach to the thermodynamic limit. To quantify 
this feature more precisely, we extrapolate our results to 
the large-TV limit. Note that we are approaching that 
limit by first taking the large-volume limit, i.e., the ther¬ 
modynamic limit is reached along the line of vanishing 
density. This is the preferred path for lattice calcula¬ 
tions in order to avoid finite-range effects. Our results 
are shown in Fig. [3j where the fits are of the form 

/(7) = / 0 + V, ( 19 ) 

and where f 0 = 1 for the energy fits and 0 for the con¬ 
tact fits. We find A = —0.60(2), B = 1.54(4) for the 
ground-state energy; A = —0.50(2), B = 1.64(4) for the 
excited-state energy; A = 0.063(3), B = 2.58(4) for the 
ground-state contact; and A = 0.041(1), B = 3.04(1) for 
the excited-state contact. In the thermodynamic limit, 
the system is known to be gapless, such that we expect 
our extrapolations of E 0 and E 1 to coincide; they do so 
only up to 7 ~ 3.0 within the extrapolation uncertain¬ 
ties. In the same figure, our answers are compared with 
those obtained directly in the thermodynamic limit in 
Ref. [Hj, which uses the Bethe ansatz. The agreement is 
not perfect but quite satisfactory, especially for 7 < 3.0. 


V. SYSTEMATIC EFFECTS 


B. Ground- and excited-state contacts 


In order to compute the contact, we proceed as with 
the energy. We use the observable of Eq. (16), which 


we can access by differentiation of In Z with respect to 
the coupling, as explained above. The interaction energy 
was thus estimated as a function of f3e F and extrapo¬ 
lated to the ground state. The extrapolation procedure 
is depicted in Fig. [5] 

From the above extrapolation, our results for the con¬ 
tact are shown in the bottom panel of Fig. [2] for a repre¬ 
sentative system. The ground-state contact values were 
taken, in a similar manner as for the ground-state energy, 
from the asymptotic behavior of the estimator for the 


In this section we elaborate on the systematic effects 
due to finite lattice size. At fixed particle number, in¬ 
creasing the lattice volume reduces the density. Thus, 
lattice-spacing effects are systematically reduced as the 
interparticle spacing grows well beyond the lattice spac¬ 
ing. The effect of increasing the lattice size is clear in 
Fig- ED where we show the energy shifted by the binding 
energy, which yields the energy of N/2 (effective) bosonic 
dimers: 


m = 3L + ^. (20 ) 

E f g E fg 7 r 

As evident in that figure, the change due to increased 
lattice size is monotonic and essentially converged at 
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Table I. Ground- and excited-state energies (E 0 and E 1 , respectively) as functions of the dimensionless coupling 7 , in units 
of E fg , extrapolated to infinite volume. The uncertainties reported are associated with the extrapolation procedure, which 
accounts for the statistical uncertainties in the Monte Carlo calculation at each finite volume. 


N 

4 

6 

8 

10 

12 

7 

E 0 

E 0 

E 1 

E 0 

E 1 

E 0 

E 1 

E 0 

E 1 

0.00 

1.5 

0 . 888 ... 

1.555... 

1.125 

1.59375 

0.96 

1.2 

1.0555... 

1.25 

0.20 

1.34(1) 

0.771(6) 

1.542(3) 

0.99(2) 

1.58(3) 

0.83(3) 

1.19(1) 

0.923(4) 

1.23(1) 

0.40 

1 . 22 ( 1 ) 

0.63(1) 

1.50(1) 

0.852(6) 

1.54(7) 

0.69(7) 

1.14(3) 

0.794(2) 

1.18(2) 

0.60 

1.07(1) 

0.51(2) 

1.43(2) 

0.711(1) 

1.5(1) 

0.55(1) 

1.08(6) 

0.649(5) 

1.10(4) 

0.80 

0.94(4) 

0.33(2) 

1.33(3) 

0.542(5) 

1.4(1) 

0.40(2) 

0.99(8) 

0.49(1) 

1 . 00 ( 6 ) 

1.00 

0.67(3) 

0.17(2) 

1.21(4) 

0.37(2) 

1 . 2 ( 2 ) 

0.23(2) 

0.9(1) 

0.31(1) 

0.86(7) 

1.20 

0.50(2) 

-0.04(2) 

1.07(5) 

0.18(1) 

1 . 1 ( 2 ) 

0.06(1) 

0.7(1) 

0.11(3) 

0.71(8) 

1.40 

0.28(3) 

-0.19(1) 

0.91(6) 

- 0 . 01 ( 2 ) 

0.9(2) 

-0.17(1) 

0 . 6 ( 2 ) 

-0.06(1) 

1.0(9) 

1.60 

0 . 10 ( 2 ) 

-0.42(1) 

0.73(7) 

-0.27(1) 

0.7(2) 

-0.37(1) 

0.4(2) 

-0.30(2) 

0.3(1) 

1.80 

- 0 . 2 ( 1 ) 

-0.65(3) 

0.52(8) 

-0.50(4) 

0.5(2) 

-0.61(3) 

0 . 2 ( 2 ) 

-0.55(2) 

0 . 1 ( 1 ) 

2.00 

-0.51(5) 

-0.88(3) 

0.30(9) 

-0.75(2) 

0.3(2) 

-0.82(2) 

-0.04(20) 

-0.79(2) 

- 0 . 1 ( 1 ) 

2.20 

-0.85(4) 

- 1 . 11 ( 2 ) 

0.06(9) 

-1.05(3) 

0.09(20) 

-1.07(1) 

-0.3(2) 

-1.03(4) 

-0.4(1) 

2.40 

-0.94(7) 

-1.41(3) 

-0.19(9) 

-1.32(1) 

- 0 . 2 ( 2 ) 

-1.37(2) 

-0.5(2) 

-1.34(2) 

-0.7(1) 

2.60 

-1.43(3) 

-1.76(4) 

-0.5(1) 

-1.62(1) 

-0.4(2) 

- 1 . 68 ( 2 ) 

- 0 . 8 ( 2 ) 

-1.65(1) 

- 1 . 0 ( 1 ) 

2.80 

-1.79(6) 

-2.08(4) 

-0.7(1) 

-1.98(3) 

-0.7(1) 

-2.01(5) 

- 1 . 1 ( 2 ) 

-1.98(1) 

-1.3(1) 

3.00 

- 2 . 2 ( 1 ) 

-2.37(7) 

- 1 . 0 ( 1 ) 

-2.35(1) 

-0.9(1) 

-2.39(2) 

-1.5(2) 

-2.33(2) 

-1.7(1) 

3.20 

-2.3(1) 

-2.83(2) 

-1.4(1) 

- 2 . 68 ( 1 ) 

- 1 . 2 ( 1 ) 

-2.73(2) 

- 1 . 8 ( 2 ) 

-2.69(3) 

- 2 . 0 ( 1 ) 

3.40 

-2.9(1) 

-3.16(3) 

-1.7(1) 

-3.07(1) 

-1.49(1) 

-3.09(4) 

- 2 . 2 ( 2 ) 

-3.16(4) 

-2.4(1) 

3.60 

-3.4(2) 

-3.55(2) 

- 2 . 0 ( 2 ) 

-3.54(3) 

- 1 . 8 ( 1 ) 

-3.56(1) 

- 2 . 6 ( 1 ) 

-3.53(4) 

- 2 . 8 ( 1 ) 

3.80 

-3.8(1) 

-4.05(2) 

-2.3(2) 

-3.95(7) 

- 2 . 2 ( 2 ) 

-4.00(4) 

-3.03(10) 

-3.96(4) 

-3.2(1) 

4.00 

-4.2(1) 

-4.51(7) 

-2.6(3) 

-4.38(5) 

-2.5(3) 

-4.49(4) 

-3.45(6) 

-4.39(7) 

-3.6(2) 


Table II. Ground- and excited-state contacts (G 0 and C \, respectively) as functions of the dimensionless coupling 7 , in units 
of Lkp, extrapolated to infinite volume. The uncertainties reported are associated with the extrapolation procedure, which 
accounts for the statistical uncertainties in the Monte Carlo calculation at each finite volume. 


N 

4 

6 

8 

10 

12 

7 

C 0 

C 0 

Cx 

C 0 

C\ 

C 0 

Ci 

C 0 

Ci 

0.00 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.20 

0.00221(7) 

0.00167(7) 

0.00023(4) 

0.00186(3) 

0.0006(4) 

0.00179(4) 

0 . 0002 ( 1 ) 

0.00185(5) 

0.0003(1) 

0.40 

0.0081(3) 

0.0073(4) 

0.0018(3) 

0.0079(2) 

0.003(2) 

0.0076(2) 

0.0014(7) 

0.00750(4) 

0 . 0020 ( 8 ) 

0.60 

0.0189(3) 

0.0169(8) 

0.007(1) 

0.0185(2) 

0.007(5) 

0.0179(5) 

0.006(3) 

0.0182(2) 

0.008(2) 

0.80 

0.034(2) 

0.034(1) 

0.018(2) 

0.0357(4) 

0.017(9) 

0.034(1) 

0.017(6) 

0.0344(9) 

0.020(3) 

1.00 

0.063(2) 

0.056(2) 

0.036(4) 

0.060(1) 

0.03(1) 

0.058(1) 

0.03(1) 

0.0594(8) 

0.040(5) 

1.20 

0.094(3) 

0.089(2) 

0.063(5) 

0.092(1) 

0.06(2) 

0.088(1) 

0.06(1) 

0.091(2) 

0.069(7) 

1.40 

0.138(3) 

0.127(1) 

0.100(7) 

0.132(2) 

0 . 10 ( 2 ) 

0.131(1) 

0 . 10 ( 2 ) 

0.132(1) 

0.109(9) 

1.60 

0.188(4) 

0.180(1) 

0.148(8) 

0.195(1) 

0.14(2) 

0.184(1) 

0.15(2) 

0.186(2) 

0.16(1) 

1.80 

0.26(2) 

0.243(3) 

0 . 21 ( 1 ) 

0.260(6) 

0 . 20 ( 1 ) 

0.249(3) 

0 . 21 ( 2 ) 

0.254(3) 

0.23(1) 

2.00 

0.35(1) 

0.317(3) 

0.28(1) 

0.342(3) 

0.264(4) 

0.326(4) 

0.29(2) 

0.335(2) 

0.31(1) 

2.20 

0.45(1) 

0.411(8) 

0.37(2) 

0.448(6) 

0.344(8) 

0.423(1) 

0.38(2) 

0.434(5) 

0.41(1) 

2.40 

0.53(1) 

0.530(9) 

0.47(2) 

0.566(1) 

0.44(2) 

0.543(2) 

0.50(2) 

0.551(3) 

0.53(2) 

2.60 

0.72(3) 

0.67(1) 

0.60(3) 

0.705(2) 

0.54(5) 

0.684(5) 

0.63(1) 

0.693(1) 

0.67(2) 

2.80 

0.90(1) 

0.83(1) 

0.73(5) 

0.879(7) 

0 . 66 ( 8 ) 

0.852(8) 

0.79(1) 

0.857(2) 

0.83(3) 

3.00 

1.07(4) 

1 . 00 ( 2 ) 

0.88(7) 

1.072(1) 

0 . 8 ( 1 ) 

1.04(1) 

0.97(2) 

1.046(3) 

1.01(4) 

3.20 

1.19(4) 

1.26(2) 

1.05(9) 
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Afj ,=100 to within the statistical uncertainty of our cal¬ 
culations. A slower rate of convergence is nevertheless 
observed as N is increased. The fact that a single uni¬ 
versal curve is achieved for each value of N indicates 
that the continuum limit is reached and is consistent with 


7 = g/n being the correct physical coupling. In order to 
extrapolate to the continuum limit, we found it sufficient 
to perform a linear fit to our data as a function of 1 /N x 
at constant N and g. 

The top panel of Fig. [4] (N= 4) also gives insight into 
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Figure 3. (Color online) Ground- and excited-state en¬ 
ergy (top) and contact (bottom) as functions of the di¬ 
mensionless coupling 7 , extrapolated to the thermodynamic 
limit. The crosses (indicated as BZ in the legend) show the 
thermodynamic-limit data from Ref. H3I. and the solid lines 
are fits (see text). 


the nature of the effective dimer-dimer interaction, which 
has been studied before with the aid of the Bethe ansatz 
and other methods (see e.g. Refs. 23 HD])- It is well 
known that, in the strong-attraction limit, the fermion 
pairs become effectively hard-core bosons 22]. As we 
see in Fig. [4] the effective boson interaction must be at 
least partially attractive in character (in that increasing 
the coupling lowers the total energy) up to 7 ~ 2 . 2 , and 
saturates or becomes mild and fully repulsive beyond that 
point. 

While the behavior of the ground-state energy with the 
lattice size is rather benign, the excited-state energies are 
affected in a more pronounced way, as shown in Fig. [5| As 
can be appreciated in that figure, the lattice effects are 
larger, which makes their extrapolation to the infinite- 
volume limit more challenging. 


VI. SUMMARY AND CONCLUSIONS 
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We have presented lattice Monte Carlo results for the 
ground- and excited-state energies as well as the ground- 
and excited-state contacts for few-fermion systems in a 
periodic one-dimensional box. Our calculations were per¬ 
formed on lattices of size N x = 20, ..., 100 and are ex¬ 
act up to statistical and systematic uncertainties, both 
of which we have addressed explicitly. Our results cover 
unpolarized systems with particle numbers in the range 
N = 4-12, and attractive coupling strengths 7 = 0, 0.2, 
..., 4.0. Although these systems can be studied directly 
with the Bethe ansatz in the continuum and thermody¬ 
namic limits, we provide here results for finite systems, 


Figure 4. (Color online) Ground-state energy of N = 4, 6 , ..., 
12 particles (top to bottom panels, respectively) in a periodic 
ID box with volumes N x = 20, 30, 50, 80, and 100. The 
binding energy of N/2 pairs, i.e. E B = Ne B /2 was added, 
where E B /E FG = 3^ 2 /- 7 T 2 and e B = 1/oq = g 2 / 4. 


spanning the few- to many-body regimes and showing ex¬ 
plicitly the approach to the thermodynamic limit for the 
energy and the contact. Our results indicate that that 
limit is approached surprisingly quickly. 

When analyzing the numerical data for excited-state 
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Figure 5. (Color online) Estimates of the ground- and excited- 
state energies of a 12-particle system in a periodic ID box with 
volumes N x = 40, 60, 80, and 100. 


quantities, we found it necessary to apply a simple yet 
powerful technique to determine the uncertainties in a 
realistic fashion. This approach consisted of performing 
power-law fits that took into account the exact results in 
the noninteracting limit. Such a constraint, together with 
the assumption that the underlying curve was smooth, 
allowed us to provide reasonable and well-defined esti¬ 
mates for the uncertainties for E 1 and C 1 , except for the 
N = 4 case, which proved numerically more challenging. 
We defer further study of that case to future work. 

Where possible, we compared our answers with pre¬ 
vious results in the thermodynamic limit (e.g., from 
Ref. 113| ) and found very good agreement. Further, our 
results serve as a benchmark for higher-dimensional stud¬ 
ies that make use of the lattice Monte Carlo method to 
study few-fermion systems at zero temperature. 

Our calculations provide a picture of the effective in¬ 
teraction between bosonic pairs in the N > 4 system. 
The nature of that interaction depends on the coupling 
strength; it lowers the energy relative to the total binding 
energy up to 7 ~ 2 . 2 , and saturates or becomes weakly 
repulsive beyond that. 

Finally, it should be stressed that the limitations we set 
in our study (namely a maximum particle number N = 
12 and lattice size N x = 100) do not reflect a breakdown 
of the method of any kind. Higher particle numbers and 
stronger couplings do require larger lattices in order to 
achieve the continuum limit, but this can also be reached 
confidently by extrapolation from smaller systems and 
lattices, as we show here. However, a better route may 
be the use of improved actions and operators (see, e.g., 
Ref. [53]) to reduce lattice spacing effects at fixed volume. 
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